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The critical behavior of a 3D Ising-like system is studied at the microscopic level of consideration. The free en- 
ergy of ordering is calculated analytically as an explicit function of temperature, an external field and the initial 
parameters of the model. Within a unified approach, both Gibbs and Helmholtz free energies are obtained and 
the dependencies of them on the external field and the order parameter, respectively, are presented graphi- 
cally The regions of stability metastability and unstability are established on the order parameter-temperature 
plane. The way of implementation of the well-known Maxwell construction is proposed at microscopic level. 
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1. Introduction 

Basic principles of phenomenological theory for second order phase transitions (PT) were formu- 
lated by Landau in late 30-40s i2j , originally to describe superconductivity. The key assumption 
of the theory is that in the vicinity of the critical point, the free energy can be expanded in a 
power series in the order parameter, the equilibrium value of which is found to obey the minimum 
condition for free energy. As is now known [3 , this expansion is well defined in the neighborhood 
of Tc (critical value of temperature) except for a narrow interval determined by the Ginzburg 
criterion [4|. In this interval, the crucial role is played by the order parameter fluctuations while 
Landau's theory is essentially a mean-field theory. Nevertheless, it is of great use as a qualitative 
tool for understanding the nature of PT. 

The above mentioned power series expansion is often called Landau's free energy. It turns 
out that similar quantity can also be derived (not constructed) in a theory considering PT at a 
microscopic level. Indeed, in the collective variables (CV) method proposed in [SJ |B] to describe 
the critical behavior of spin systems, such a quantity appears in a natural way. In particular, a 
microscopic analogue of Landau free energy was calculated in |7j, where the Ising-like system was 
considered in zero external field. 

The purpose of the present paper is to extend the results obtained earlier by the CV method for 
physical characteristics of 3D Ising-like model to the presence of an external field and look at what 
is going on when the field changes sign. As a result, we obtain the Gibbs free energy, Helmholtz free 
energy, and the order parameter and establish the regions of stability, metastability, and unstability 
for the system under investigation. 



In our research we consider the system of N Ising spins placed on sites of a simple cubic lattice 
of the spacing c. The Hamiltonian of such a system in an external field is well known 



2. Method 




(2.1) 
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Here, is a short-range interaction potential between spins located at the i-th and j-th sites, 

the spin variables cr.; take on values ±1, H is the external field. We do not restrict the summation 
in ( |2.1[ ) to the nearest neighbors. To this end, <i>(r) = const x cxp(— r/6) can be chosen as the 
interaction potential with an efi^ective range b. 

It is also known that this problem has not yet been solved. Universal critical characteristics of 
three-dimensional (3D) systems are successfully described in a variety of approaches. Among them 
it is worth mentioning the field theoretical and renormalization group (RG) methods [5^ , which also 
enable one to calculate nonuniversal characteristics (critical amplitudes, critical temperatures), 
but without taking into account the dependency on initial parameters of the Hamiltonian. An 
alternative method for theoretical investigation of the model under consideration is the method of 
CV [5_, which has an advantage of being a successive microscopic approach, which, in turn, makes 
it possible for thermodynamic functions and physical characteristics of the system to be explicitly 
obtained as functions of temperature, the field, and microscopic parameters of the model. 

In the framework of the "p'*-model" approximation the functional representation for the parti- 
tion function in terms of CV pk is as follows: 

Z = Zq l{dp)^°exp aiy^po - ^ ^ d(A;)pkp_k - 



4! 



^ Pki •■•/3k4'5ki + ...+k4 



(2.2) 



kiSBo 

Here, the quantity d{k) contains the Fourier transform of the interaction potential 

d(/c) =a2+/3$(0)$-/3$(A:). (2.3) 

The explicit expressions for a„ can be found in j9j (equation (1.25)) , iVo ~ N/s^ where the 
quantity sq defines the region of validity for the parabolic approximation of Fourier transform of 



the interaction potential (for details, see [TU]). The summation in (2.2) is performed over wave 
vectors of the first Brillouin zone corresponding to a reciprocal effective lattice with the lattice 
constant csq. 



In the work jTTj the method for calculation of the partition function (2.2 1 of the model near 
the second order PT point was generalized to the case of the presence of the fixed external field of 
an arbitrary magnitude. Unlike the work [12], here no assumption concerning either the strength 
or the weakness of the applied field is made and, therefore, no perturbation series are used. The 
calculation procedure is based on Kadanoff's idea of constructing effective spin blocks [T3]. We 
divide the phase space of the CV pk into layers, the RG parameter being s, and average the 



Fourier transform of the potential in each n-th layer. The first step in calculating equation (2.2) is 



to integrate over pk with k € [fcmax, ^max/s], ^max = tt/cs q. T he result of integration proves to be 



represented in the form of the initial expression, equation (|2.2|), with renormalized coefficients a'"^\ 



di{k), and a^^\ If we perform a step-by-step integration of the partition function over rip layers, 
we arrive at 



Z = Zo[Q{d)f°( YlQnjZ^GR 



(2.4) 



(n) 

The partial partition functions Qn of the n-ih level are characterized by a set of coefficients , 
dn{0), af'\ for which the recurrence relations (RR) hold (the work [9] is devoted to this problem). 
The quantity Up is called the exit point from the critical regime of the order parameter fluctuations. 
It defines the number of iterations at which the system is still in the scaling region. 

In [T^, two regions in (/i, r)-plane were distinguished, which are of weak and of strong field. The 
calculations in the regions were performed using different forms of exit points. These quantities 
were subject to certain conditions which we do not discuss in the present work. If the field was 
considered to be strong, the formula Up = rih = — In /i/ In _Bi — 1 was used, while in the case 
of a weak field another expression Up — = — ln-r/lni?2 — 1 was taken. Hence, the question 



43002-2 



Free energy of an Ising-like model 



immediately arises what should be done at intermediate values of a field. The problem is solved 
by constructing a new expression for the exit point, proposed in [llj : 



In U 



21n£;i 



(2.5) 



where some temperature fields hi'^'' = \t\p° and hi ' = Irij^" are introduced, po = ln£'i/lni?2 
is the so-called "gap exponent" [M], the signs "+" and "-" are related to T > Tc and T < Tc 
respectively, and 

h^l3n, r^^T, T, = -ErT, r=^^. (2.6) 

"0 Jo J c 

The quantities Ei and E2 are the eigenvalues of the matrix of the RG transformation linearized 
near the fixed point of RR. 



In the limiting cases, as /i or as r — > 0, equation (2.5 1 takes on the form of to^ or n/i, 
respectively, and therefore can be applied to the case of an arbitrary field. Nevertheless, it is worth 
mentioning that the inequality h ^ he defines a strong field, and h <^ h^ defines a weak field. 
Note also that there is an arbitrariness in choosing rip, which is discussed in [TUIITT] in more detail. 
What is important is that for n > Up one should be able to integrate the partition function with 
Gaussian measure density. 



Table 1. Numerical values of the parameters used in present calculations. 



b/c /?c$(0) so 


/o 


ho 


Cfci 


00 




no 


0.3 1.6411 2.0 


0.5 


0.760 


1.1762 


0.5938 


0.105 


0.5 



We will adhere to the approach adopted in [TS], where the coefficients c^^ , fo, ho are also defined 
(see equations (4.4), (4.17), (4.23), and (4.49) in [H]), and take Ei = 24.551, E2 = 8.308. The 
information on the physical meaning of coefficient no can be found in |16| . Numerical values of all 
coefficients used to obtain the final graphical results are presented in table [l] It is worth stressing 
that having fixed the value of RG parameter s = s* {s* = 3.5977 in the present calculation), there 
remains only one initial parameter, which is the ratio of the effective range of interaction b to the 
lattice constant c. All the other coefficients can be expressed via b/c [TUl [T7] . 



3. Results and discussion 

Based on the above mentioned method, in the work [TT] as well as in jTH] free energy of 3D 
Ising-like model was calculated as a logarithm of the partition function (2.4) multiplied by —kT 
and presented in the form of several contributions 

F{T,h)^Fa+F^^^+Fl,^\ (3.1) 

Here, the term Fa is the analytical part of free energy and does not affect the critical behavior of 
the system. The last two terms in r.h.s. of (3.1 1 have non-analytical dependence on temperature 
and on the external field. The explicit expressions for -F'i^'' and i^o^"* are 

(r, h) = -kTN^i^) {h^ + h^y') ^ , (3.2) 

and 

Fl,^\T,h) = -kTNEo{<y±). (3.3) 
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The quantity 7!*^ from (3.2 1 includes contributions from the critical regime of the order parameter 
fluctuations and from the limiting Gaussian regime (LGR)(for details, see [lliilBj ). The mentioned 
critical regime is characterized by the RG symmetry. 

The quantity Eo{(J±) is the contribution from the collective variable pQ. As is known from the 
theory of CV [5 , the mean value of pk=o is connected with the order parameter and consequently 
the quantity Fq is the free energy of ordering [7]. If there is no external field, Fq is analogous to 
Landau's free energy in phenomenological theory of second order PT [7J. Since we work in the 
ensemble of N particles, with field h and temperature r being independent variables, this analogy 
is now not direct. Following Stanley [14 , and based on convincing arguments of the work [18] , we 
associate the partition function of our model with Gibbs free energy. An appropriate thermody- 
namic potential is Gibbs free energy provided the magnetic field and temperature are independent 
("natural") variables. If the role of independent variables is played by an order parameter (here 
magnetization per spin) and temperature, then the thermodynamic potential associated with the 
partition function is Helmholtz free energy. Hence, Landau's free energy is rather Helmholtz than 



Gibbs one. The expressions (3.1 1-(3.3 1 in turn present Gibbs free energy for the considered model, 
although this was not specified in earlier works, particularly in [IH I16j. Therefore, in order to ob- 
tain a microscopic analogue of Landau's energy one should perform the Legendre transformation. 
Let ^0 denote this analogue. Then, 

AoiT,M)=Fo + Mon, (3.4) 
where Mq is the magnetization of the system under consideration such that 

1 dFo dEo 

The expression for i?o(o'±), which was found in [TT] (see equation (4.4) there) for T > and 
in |T5] (see equation (4.3) there) for T < Tc, reads 

1 3 3(np-|-2) 

Eoi<y±) = ha± - -d„^+2(0)4 - ai"-+'Vi . (3.6) 

The quantity a± was found from condition 

dEo{a± 
da± 

in the form 



(3.7) 



a± (3.8) 
This results in the cubic equation for ctq*'' 

al+p<Jo + q^O (3.9) 



with coefficients 



= -6.0^/^^^/^^^^^ ^-T7^, (3.10) 



where the denotations r„ and m„ are connected with (i„(0) and 04 through (i„(0) — s 
s^^Un, and according to [ITJ[T2] are expressed as 

rn,+2 = /3c$(0)/o (-1 ± f i?2^+' 



Un,+2 = WMO)Y^o[l±fE'^''^'^f\. (3.11) 
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Figure 1. (Color online) Solutions of cubic equation (3.91 (a 
(b) as functions of field for r — ±10~^ 



h = 10' 



real solutions at r 0; marked by "01", "02", and 
respect, we will refer to the solutions as a, 



(+) 



as functions of temperature for 
The curves marked by "0" correspond to the 
03" denote different solution at r < 0. In this 
1,2,3). 



and cr^. > {i = 



Note, that the coefficients p, q should be marked by the superscript ± (as should be done for Up 
from (2.5 1), but where it does not cause confusion, we drop this superscript out. One should just 
remember about a difference in temperature scales between the cases of T > Tc and of T < Tc. 



Equation (3.9 1 can be solved by Cardano's method. The solutions are presented graphically 



in figure [T] From the plot we see that for a given value of field there exists some temperature 
To < such that equation (3.9 1 has three real solutions for r ^ tq and one real solution for r > tq. 
This quantity, tq, corresponds to the condition Q = where Q is the discriminant of the cubic 



equation (3.9): 



(9/2)' 



(3.12) 



with q and p from (3.10 I. When h ^ 0, then Tq = 0. If h ^ 0, then Tq is determined by setting r.h.s. 
of (3.12) equal to zero. In [16 , tq was found numerically and the plot of tq versus h was drawn. 
As we will see in what follows, the solution tq = To(h) to the equation Q = defines a curve in 
"order parameter-temperature" plane which can be identified with the spinodal of a fluid. 



Taking into account (2.5), and (3.8), the quantity -Eo(f) takes on 



where we have introduced the following notation 



(3.13) 



(±) 



3^(±)- 



(3.14) 



In the previous works |10l [TTl [T51 [TB] the system was considered in the external field h ^ 0. 
Analytical results for (Gibbs) free energy [TT| ITS] , the order parameter |10[ [TB] , and the suscepti- 
bility [1^ were obtained. The purpose of this paper is to extend the method to the region h < 
and to take into account all solutions to the equation (3.9). We look at Gibbs free energy Fg and 
at the contribution from it to the order parameter Mq of equation ( |3.5[ ). Thereupon, Hclmholtz 
free energy is calculated. 

The main results of this paper are presented in figures [2] [3] and|4] In figure |2] we can see: (a) 
Gibbs free energy Fq as a function of the external field in low temperature region (t — —0.001); 
(b) a plot of Fq versus Mq. Each of the points a, b, c, . . . represents a certain specific state of 
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Figure 2. (Color online) A set of 4 pictures is placed for the sake of clarity. All are drawn for 
T = —0.001. The top left is the field dependence of Gibbs free energy Fq. The top right shows Fq 
against A4o- The bottom left presents the solutions of the cubic equation | |3.9[ ). The bottom right 
is the equation of state for the system under consideration. The values of energies are normalized 
by dividing by kTN. The points a to / represent some particular states of the system in different 
coordinates. See the text for details. 




Figure 3. (Color online) The plot shows temperature dependence of magnetization calculated 
at different solutions to equation ( |3.9[ ). Thick line is magnetization in zero external field and 
corresponds to the coexistence curve (binodal). In negative external field, the solution croi gives 
rise to Curves 1-5 on the magnetization-temperature plane. In positive external field, ctos gives 
rise to Curves 6-10. Curves 11-20, which correspond to ao2, can be interpreted as non-physical. 
Thick dashed curve corresponds to the saturation curve (spinodal). In order to recover that, one 
should compute magnetization corresponding to (tqi at h < and tq = ro(|/i|) (upper branch) 
and corresponding to aoa at /i > and tq = To{h) (lower branch). The magnitudes of field, at 
which the results are drawn, are \h\ = 0, 10"^, 10"®, 5 • 10"®, 10"^, 2 • 10"^ 
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Figure 4. (Color online) Helmholtz free energy Aq as a function of order parameter for three 
different temperatures T < Tc {r = -0.001), T = Tc {t = 0), and T > Tc {r = 0.001). The 
value is normalized by dividing by kTN . 



the system in different coordinates. Especially, points a, c, d, and / correspond to such a relation 
between the field h and temperature r that the equality r = to(/i) holds. Thus, in each of these 
we have Q = Q. Points b and e are associated with first order phase transition and correspond 
to /i = and r = —0.001. The isotherm in figure |2] (b) was drawn with the help of parametric 
representation {Fq ~ FQ{h), AIq = Mo(/i)}. 

Recall that Gibbs free energy is a concave function of field [14J. Therefore, from figure l2| (a) we 

(— ) . . . . 

can conclude that the solution cTgj gives rise to a non-physical result. It corresponds to the region 

of unstability {c—d part of the curve in figure 2 (a), (b)). This is not surprising if we note that 



02 

maximizes Fg. Solutions ctq]^ •* at /i < and CTpg^ at /i > correspond to the region of metastable 
states (d — e and b — c parts of the curve respectively) . In the case of liquid-gas PT it would be 
the region between the binodal and the spinodal curves. Hence, we regard these results important 
from the perspective of using this method for description of the critical behavior in simple fiuids. 

Figure [2] (c) to some extent repeats figure [l] (b) and is placed here for convenience of repre- 
sentation. The quantity ctq is related to the scaling function of magnetization via formula (4.8) in 
work [16]. Finally, in figure [2] (d) the equation of state computed by means of ( |3.5| is presented at 
T = -0.001. 

A more detailed representation of the equation of state is shown in figure [3] The derivatives of 
Gibbs free energy with respect to field are presented as functions of temperature, different solutions 
to equation (3.9 1 being taken into account. As we can see, such an accounting allows us to establish 
spinodal and binodal curves of the model and, in this respect, to find the stable, metastable, and 
unstable regions on "order parameter-temperature" plane. 

The order parameter dependence of free energy Ag is presented in figure |4] for three different 
temperatures T > Tc, T — Tc and T < Tc- Here again a, b, c, . . . denote the states as in figure|2] We 
can see that Helmholtz free energy has two minima below Tc- The equilibrium values for the order 
parameter are defined by these minima. Knowing their positions at different temperatures, we can 
recover the coexistence curve. Points b and e coincide with the minima at t = —0.001. Note that 
based on the principle of minimum Gibbs free energy and from figure [2] (a), we can conclude that 
with the field changing from h < to h > the system will tend to move along a — b — e — f part 
of the curve. Applied to Helmholtz free energy, it means that the system will tend to jump from 
state b to state e. Therefore, one may supplement figure |4] with a double-tangent construction. As 
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a consequence, on Mq — h plane we have a horizontal segment for Mq at h = (see figure |2] (d)), 
which corresponds to a Maxwell construction derived at microscopic level. On the other hand, 
based on the maxima of Gibbs free energy with respect to Mq (see figure [2] (b)), we are able to 
construct a saturation curve. Let us just remember that dependence Fq on the order parameter is 
formal due to AIq being not "natural" variable for Fq. However, the same results for the binodal 
and the spinodal can immediately be obtained if one appropriately chooses the solutions to (3.9). 

The forms of isotherms in figure |4] are similar to those in Landau's theory [3 . However, the 
peculiar feature is that in the present work the free energy is an explicit function of temperature, 
external field and microscopic parameters (in the present case a ratio b/c) of the model, with a 
non-analytical dependence on its arguments. This gives rise to a non-classical critical behavior of 
the system, i.e., to the critical exponents taking on non-classical values. These values for the most 
important exponents are reported in table [2j There are also collected the values obtainable by 
CV method in higher, approximation as reported in Conclusions of the work f7^. Note that in 
our calculation we have neglected the critical exponent rj responsible for the behavior of the pair 
correlation function at T — Tc. Accounting for corrections to scaling is beyond the scope of the 
present paper as well. In table [2] the classic values for the critical exponents are presented as well 
as the results of other authors. We stick to the following notation. The temperature behavior of 
magnetization is governed by /3, of the heat capacity by a, of susceptibility by 7, of the correlation 
length by ly. The field behavior of magnetization is governed by S, M ^ \h\^^^ . We have calculated 
some of the exponents with the help of scaling laws. Based on earlier results [15_, we are also able 
to estimate the critical exponents (p and ip that describe the field behavior of the heat capacity, 
C ^ and of the entropy, S ^ \h\^ , respectively. Their numerical values are (p = 0.122 and 

ip = 0.539, which differ from tp — and -0 = 2/3 in mean-field theory |14j . 



TablS 2. Numerical values of the critical exponents. CV, collective variables method; MF, mean- 
field values; HT, high-temperature expansion results; FT, field theory approach. 





CV, 
p*-model 


CV, 
p^-model [7| 


MF 




FT lU 


a 


0.185 


0.088 





0.110(1) 


0.109(4) 


/? 


0.302 


0.319 


1/2 


0.3263(4) 


0.3258(14) 


7 


1.210 


1.275 


1 


1.2373(2) 


1.2396(13) 


V 


0.605 


0.637 


1/2 


0.6301(2) 


0.6304(13) 


S 


5 


5 


3 


4.792 


4.805 



It should also be mentioned that the scaling properties for physical characteristics of the con- 
sidered model are discussed in earlier works. In section 3 of the work [TD] the scaling functions for 
(Gibbs) free energy, for the order parameter, and for susceptibility were calculated. 



4. Conclusions 

In this work, an analytical expression for Gibbs free energy of Ising-like system is obtained and 
investigated near the second order phase transition. The emphasis is made on the presence of the 
external field and on the situation where the field changes its direction. Some well-defined concepts 
of phenomenological theory of phase transitions - e.g. Landau's energy, the Maxwell construction, 
the double-tangent construction - are derived and reexamined at microscopic level. We hope that 
the approach described and the results presented should be useful in further investigations of 
critical behavior in 3D Ising-like systems with analytical methods; in particular, in order to obtain 
the coexistence curve and the spinodal curve for a system at least in the vicinity of critical point. 
There exists a general belief |14[ 1^ [22] that simple fluids belong to the Ising universality class. 
Hence, we regard these results important from the perspective of using this method for description 
of the critical behavior in simple fluids. 
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Bi/ibHa eHeprifl ri66ca Ta ei/ibHa eHepria Fe/ibMro/ibi^a fl/iJi 
Tpi/iBMMipHoV i3iHronoAi6HoV MOA^-ni 

M.n. Kos/ioBCbKMM, P.B. PoManiK 

iHCTHTyT 4>i3i/iKi/i KOH/^eHCOBaHi/ix CMCreM HAH VKpaiHi/i, Byn. I. CBeHL^iLfbKoro, 1, 79011 JlbsiB, VKpaiHa 

B flaHiw po5oTi Ha MiKpocKoninHOMy piBHi poBr/isqfly BWBHaeTbca Kpi/iTWHHa noBefliHKa rpuBi/iMipHoT isiHro- 
nofliSHoT cwcTeMW. AHa/iixwHHO o6HHC/iK)eTbCii Bi/ibHa eHepria BnopaflKyBaHHa jik 4>yHKLi,i!R TeMnepaxypi/i, 
BOBHiujHboro HOJiiR i noHaxKOBMX napaMexpiB MOfle/ii. B Me>Kax sanponoHOBaHoro niflxofly oxpuMaHO Bnpasw 
fl/ia Bi/ibHoT eHepriV ri56ca i Fe/ibMro/ibLia, Tx 3a/ie>KHOcxi Bifl nonsf i napaMexpa nopaflKy npwBefleHi xaKO>K 
rpacjjiHHO. Ha n/iomwHi napaMexp nopjiflKy-xeMnepaxypa BHaCifleHi o5;iacxi cxa5i;ibHocxi, Mexacxa6i;ibHOCxi 
xa Hecxa6mbHocxi. SanponoHOBaHO cnoci6 peajii3aL4iT npaBi/iyia MaKcseji/ia Ha MiKpoCKOninnOMy piBhi. 

K/iioHOBi c/iOBa: MOfleJib Isinra, Kpi^m^Ha noBeft'iHKa, soBHiiiiHe nojie 
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